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Abstract 

We show how one can actually take advantage of the strongly non- 
Gaussian nature of the fluctuations of financial assets to simplify the 
calculation of the Value-at-Risk of complex non linear portfolios. The 
resulting equations are not hard to solve numerically, and should allow 
fast VaR and AVaR estimates of large portfolios, where by construc- 
tion the influence of rare events is taken into account reliably. Our 
method can be seen as a correctly probabilized 'scenario' calculation 
(or 'stress-testing'). 



1 Introduction 

A very important issue for the control of risk of complex portfolios, which 
involves many non linear assets such as options, interest rate derivatives, 
etc. is to be able to estimate reliably its Value-at-Risk, or equivalently the 
probability of large downward moves, deep in the tails of the probability 
distributions ]l|, ||. This is a difficult problem, since both the non-Gaussian 
nature of the fluctuations of the underlying assets and the non-linear de- 
pendence of the price of the derivatives must be dealt with. A solution to 
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cope with non-linearity is to use Monte-Carlo simulations based on Gaussian 
multivariate statistics for the time evolutions of all the assets underlying the 
portfolio. This solution is however time consuming (specially to obtain good 
statistics in the tails of the distribution) and cannot be used for fast VaR, 
or AVaR calculations, which are important for real time estimates of the 
influence of a particular trade on the global exposure of a portfolio. More 
importantly, this method is not reliable because of the strongly non-Gaussian 
nature of the extreme moves. 'Fat tails' effects are well known and lead to 
a significant increase of the VaR estimate, even in the simplest case of a 
linear portfolio, for example containing stocks only. These fat tails can be 
further amplified by the non linear nature of the relation between derivative 
products and the underlying assets, thereby leading to very large differences 
between a Gaussian VaR estimate and reality. 

The aim of this paper is to introduce a method, called the 'optimal fluc- 
tuation method' in the physics hterature 0. This method is well suited 
to estimate large risks in the case where the fluctuations of the 'explicative 
factors' are strongly non-Gaussian (a more precise statement will be made 
below). An approximate formula can be obtained for the Value-at-Risk of a 
general non linear portfolio. This formula can easily be implemented numeri- 
cally. The basic idea is to identify the 'most dangerous' market configuration 
for a given portfolio. In the case of strongly non Gaussian fluctuations, the 
largest moves of the portfolio correspond to a large change in one explicative 
factor, accompanied by the simultaneous 'typical' changes of all the others. 
In a Gaussian world, on the opposite, the large moves of a portfolio cor- 
respond to a 'conspiracy', where all factors coherently change by a small 
amount. 

2 Fat-tailed explicative factors 

Let us assume that the variations of the value of the portfolio can be written 
as a function (i/(ei, 62, cm) of a set of M independent random variables 
Ca, a = 1, M, called 'explicative factors'. These factors can be determined 
by a classical Principal Component Analysis, where the correlation matrix of 
all assets' increments is diagonalized. However, it might be more useful for 
our purpose to consider other definitions of the correlation, more suited to 
tail events (see, e.g. p|). On short time scales, most relevant for VaR esti- 
mates, all trend effects are negligible, and we shall therefore set (ca) = and 
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{GaGb) = ^a,60'a! whciG (...) dcHotes the averagc over the relevant probabihty 
distribution. The sensitivity of the portfoho to these 'exphcative factors' can 
be measured as the derivatives of the value of the portfolio with respect to 
the Ca- We shall therefore introduce the different A's and F's as: 

df d^f 

Aa = - — Fab = — . (1) 

Ota deader 

We are interested in the probabihty for a large fluctuation df* of the port- 
folio. We will first surmise that this is due to a particularly large fluctuation 
of one explicative factor - say a = 1 - that we will call the 'dominant' factor. 
(The generalization to several factors will be discussed below). Note that 
this is not always true, and depends on the statistics of the fluctuations of 
the Ca- A condition for this assumption to be true will be discussed below, 
and requires in particular that the tail of the dominant factor should not 
decrease faster than an exponential. Fortunately, this is a good assumption 
in financial markets, but would be completely wrong for Gaussian statistics. 
Note also that the dominant factor depends a priori, via the A's, on the 
portfolios composition. 



3 The dominant factor approximation 

The aim is to compute the Value-at-Risk of a certain portfolio. This is 
defined as the value df* such that the probabihty V>{df*) (defined as the 
cumulative probability that the variation of / exceeds df*) is equal to a 
certain probability p - say 1% for a 99% confidence VaR. Our assumption 
about the existence of a dominant factor means that these events correspond 
to a market configuration where the fluctuation Sci is large, while all other 
factors are relatively small. Therefore, the large variations of the portfoho 
can be approximated as: 

M ^ M 

df{ei, 62, Cm) = df{ei) + ^ A„e„ + ^ X) ^a,beaeb, (2) 

a=2 ^ a,6=2 

where df{ei) is a shorthand notation for df{ei, 0, ....0), and all the derivatives 
are calculated at the point (ei, 0, ....0). Now, we use the fact that: 

. M 

'P>{df*) = Wdca P(ei, 62, eM)e(d/(ei, 62, cm) - df*), (3) 

a=l 
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where B(x > 0) = 1 and G(x < 0) = is the Heaviside function. We now 
expand the B function to second order, and perform the integration over the 
Ca's (a > 1), to finally obtain (see for details): 

^>(rfr) =^>(eD + E^P(eD-E^ (^P'(eD + ^P(enj , (4) 

where e\ is such that df{el) = df*, and is computed for ei = e^, ea>i = 
and where -P(ei) is the marginal probability distribution of the first factor. 
This probability density can be estimated empirically, and fitted to one of 
the possible distribution known to describe well financial data, such as a 
Truncated Levy, Hyperbolic, or Student distribution |0. 

In order to find the Value- at-Risk df*, one should thus solve (H) for el 
with V^{df*) = p, and then compute (i/(e^, 0, 0). Note that the equation 
is not trivial since the Greeks must be estimated at the solution point In 
practice, the dominant factor can be found by trial and error, by computing 
e* for all M explicative factors, and choosing the one that leads to the largest 
VaR. 

4 Discussion of the result 

Let us discuss the general result (^) in the simple case of a linear portfolio 
of assets, such that no convexity is present: the A^'s are constant and the 
Ta,aS are all zero. The equation then takes the following simpler form: 

M ^2 2 

^>(ei)-E^nen=p- (5) 

a=2 

Naively, one could have thought that in the dominant factor approximation, 
the value of would be the Value-at-Risk value of ei for the probability p, 
defined as: 

VyieiyaR) = p. (6) 

However, the above equation shows that there is a correction term propor- 
tional to -P'(e*). Since the latter quantity is negative, one sees that e* is 
actually larger than eiyaR, and therefore df* > df{eiyaR). This reflects the 
effect of all other factors, which tend to increase the Value-at-Risk of the 
portfolio. 
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The result obtained above relies on a second order expansion; when are 
higher order corrections negligible? It is easy to see that higher order terms 
involve higher order derivatives of -P(ei). A condition for these terms to be 
negligible in the limit p ^ 0, or e* — oo, is that the successive derivatives of 
P(ei) become smaller and smaller. This is true provided that P{ei) decays 
more slowly than exponentially, for example as a power-law. In this case, 
which corresponds to financial reality |@, |, |^ , each term in the expansion is 
a factor 1/e* smaller than the previous one, which indeed becomes negligible 
in the limit p —>■ 0, oo. On the contrary, when -P(ei) decays faster 

than exponentially (for example in the Gaussian case), then the expansion 
proposed above completely looses its meaning, since higher and higher correc- 
tions become dominant when p — >■ 0. This is indeed expected: in a Gaussian 
world, a large event results from the accidental superposition of many small 
events; the idea of expanding around one large even is therefore not adapted. 
In a power-law world, large events are associated to one single large fluctua- 
tion which dominates over all the others, and the above method is congenial 
for fast and precise VaR estimates. The limiting case where -P(ei) decays as 
an exponential is interesting, since it is often a good approximation for the 
tail of the fluctuations of financial assets. Taking -P(ei) ~ ai exp — a;i|ei|, 
one finds that el is the solution of: 
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ai 




(7) 



Since one has af oc af^, the correction term is small provided that the 
variance of the portfolio generated by the dominant factor is much larger 
than the sum of the variance of all other factors. 

Coming back to equation (^), one expects that if the dominant factor is 
correctly identified, and if the distribution is such that the above expansion 
makes sense, an approximate solution is given by e* = eiyaB. + with: 

where now all the Greeks at estimated at ei vaR- 
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5 Numerical tests and improved approximation schemes 

We have tested numerically the above idea by calculating the VaR of a simple 
portfolios which depend on four independent factors Ci, ...,64, that we chose 
to be independent Student variables of unit variance with a tail exponent 
fjL — A. We have considered two different portfoHos, 'hnear' (L), such that 
its variations are given by d£ = ei + + |e3 + ^64, and 'quadratic' (Q), 
such that its variations are given by dQ = dC + {dCy. We have determined 
the 99%, 99.5% and 99.9% VaR of these portfolios both using a Monte-Carlo 
(mc) scheme and the above approach. The results are summarized in Table 
1, where we give the simple estimate based on the VaR on the most dangerous 
factor eiyaR: (called Th. 0) and the estimate which includes the contribution 
of the other factors e* (called Th. 1). For both L and Q, one sees that the 
latter estimate represents a significant improvement over the naive eiyaR 
estimate. One the other hand, our calculation still underestimates the MC 
result, in particular for the Q portfolio. This can be traced back to the 
fact that there are actually other different dangerous market configurations 
which contribute to the VaR for this particular choice of parameters. Our 
formalism can however easily be adapted to the case where two (or more) 
dangerous configurations need to be considered. The general equations read: 



where a — 1,...,K are the K different dangerous factors. The e*, and there- 
fore df*, are determined by the following K conditions: 



dfiel) = dfie;) = ... = drie*^) V>,+V>2 + ... + V>k = P- (10) 



We have computed the VaR for the different portfolios in the K = 2 approxi- 
mation. The results are reported in Table 1 under the name 'Th. 2'; one sees 
that the agreement with the MC simulations is improved, and is actually ex- 
cellent in the linear case, where both the configurations ei positive and large 
and 62 positive and large contribute. In the quadratic case, the two most 
dangerous configurations are ei positive and large and ei negative and large. 
In order to get perfect agreement with the Monte-Carlo result, one should 
extend the calculation to i^T = 3, taking into account the configuration where 
62 positive and large. 




(9) 
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P = 


1% 


p = 0.5% 


p = 0.1% 




C 


Q 


C 


Q 


C 


Q 


MC 


2.93 


13.3 


3.53 


18.7 


5.30 


40.6 


Th. 


2.65 


9.7 


3.26 


13.9 


5.07 


30.8 


Th. 1 


2.83 


10.9 


3.42 


15.1 


5.20 


32.2 


Th. 2 


2.93 


12.1 


3.52 


17.2 


5.30 


38.6 



Table 1: 

6 Conclusion 

In summary, we have shown how one can actually take advantage of the 
strongly non-Gaussian nature of the fluctuations of financial assets to sim- 
plify the analytic calculation of the Value- at- Risk of complicated non linear 
portfolios. The resulting equations are not hard to solve numerically, and 
should allow fast VaR and AVaR estimates of large portfolios, where by con- 
struction the influence of rare events is taken into account reliably. This 
method is a short-cut to Monte-Carlo methods, in the sense that we directly 
identify the events which are most relevant for extreme risks, without having 
to 'wait' for them to appear in a Monte-Carlo sampling. Interestingly, the 
calculation allows one to visualize directly the 'dangerous' market configura- 
tions which corresponds to these extreme risks, by re-expressing all the real 
assets variations in terms of the dominant factors. In this sense, our method 
corresponds to a correctly probabilized 'scenario' calculation. 
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